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Abstract 

Mutation of Class bottom boat, the Drosophila homologue of the bone morphogenetic protein or growth/differentiation 
factor (BMP/GDF) family of genes in vertebrates, has been shown to disrupt development of neuromuscular junctions (NIVIJ). 
Here we tested whether this same conclusion can be broadened to vertebrate BMP/GDF genes. This analysis was also 
extended to consider whether such genes are required for NMJ maintenance in post-larval stages, as this would argue that 
BMP genes are viable candidates for analysis in progressive neuromuscular disease. Zebrafish mutants harboring 
homozygous null mutations in the BMP-family gene gdfiSo were raised to adulthood and assessed for neuromuscular 
deficits. Fish lacking gdfSa exhibited decreased endurance {~50%, p = 0.005) compared to wild type, and this deficit 
progressively worsened with age. These fish also presented with significantly disrupted NMJ morphology {p = 0.009), and a 
lower abundance of spinal motor neurons {~50%, p<0.001) compared to wild type. Noting the similarity of these 
symptoms to those of Amyotrophic Lateral Sclerosis (ALS) model mice and fish, we asked if mutations in gclf6a would 
enhance the phenotypes observed in the latter, i.e. in zebrafish over-expressing mutant Superoxide Dismutase 1 (SODl). 
Amongst younger adult fish only bigenic fish harboring both the SODl transgene and gdf6a mutations, but not siblings 
with other combinations of these gene modifications, displayed significantly reduced endurance (75%, p<0.05) and 
strength/power (75%, p<0.05), as well as disrupted NMJ morphology (p<0.001) compared to wild type siblings. Bigenic fish 
also had lower survival rates compared to other genotypes. Thus conclusions regarding a role for BMP ligands in effecting 
NMJ can be extended to vertebrates, supporting conservation of mechanisms relevant to neuromuscular degenerative 
diseases. These conclusions synergize with past findings to argue for further analysis of CDF6 and other BMP genes as 
modifier loci, potentially affecting susceptibility to ALS and perhaps a broader suite of neurodegenerative diseases. 
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Introduction 

Development, growth and stabilization of neuromu.scular 
junction.s (NMJ) in larval Drosophila require the bone morphoge- 
netic protein (BMP) gene glass bottom boat (gbb) [1]. Here we tested 
whether this same conclusion can be broadened to vertebrates, 
and extended this analysis to consider whether such genes are 
required for vertebrate NMJ maintenance in post-larval stages. 
Gbb is the Drosophila homolog of a family of vertebrate genes 
including the BMP and growth/differentiation factor ligands 
(BMP/ GDF Family), which is itself a sub-family of the transform- 
ing growth factor P (TGPP) genes [1-3]. Extending their role to 
encompass vertebrate NMJs would embolden speculation that 
mutations in BMP/ GDF genes can sensitize patients to progressive 
late-onset neuromuscular disease. This affirmation would support 
the contended relevance of several intriguing cellular and 
molecular mechanisms of neuromuscular degeneration, derived 
from Drosophila studies (see Discussion), that are hypothesized to 



impinge upon the development of treatments or diagnostics of 
neuromuscular disease. 

Indeed the synthesis of disparate literature by several authors 
has recendy suggested that mutations in BMP/GDFs are good 
candidates for sensitizing patients to amyotrophic lateral sclerosis 
(ALS), if not representing causal instigators of disease etiology 
[4,5] . This suggestion was based not only on the aforementioned 
role of Gbb [1,2], but also upon the requirement for the proteins 
constituting Gbb receptors in NMJ development [6-8]. BMP/ 
GDFs have also been implicated in ALS progression via a separate 
line of inquiry, in that a Drosophila model of familial ALS8 
demonstrates disrupted BMP signaling at their NMJs; Thus 
mutations in VAPB (vesicle-associated membrane protein B) cause 
ALS8 [9,10] and altering VapB disrupts BMP signahng at the 
Drosophila NMJ [11]. A role for BMP/GDFs in other neuromus- 
cular diseases has also been proposed, including Spinal Muscular 
Atrophy, Hereditary Spastic Paraplegias, Multiple Sclerosis and 
Huntington's Disease [4]. 
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BMPs are most widely recognized for their fundamental roles in 
development across vertebrates, including patterning the dorso- 
ventral axis of the body, CNS and retina [12-17]. BMP/GDF 
proteins heterodimerize or homodimerize to signal through BMP 
receptors, and canonicaUy through phosphorylatrng SMAD 
proteins, though several other signaling cascades can be important 
in many instances [18]. 

ALS is a progressive neuromuscular disease caused by 
motoneuron loss, though the etiology of motoneuron death is 
unknown. Candidate causes include glutamate excitotoxicity, 
oxidative stress, and RNA processing defects. Symptoms in 
patients, recapitulated in animal models overexpressing mutant 
SODl, include progressive muscle weakness and decreased 
endurance, altered gait, motor neuron death and progressive 
paralysis. Genetics of familial ALS (fALS) include lesions in 
TARDBP, FUS/TLS, C9orf72 and SODl. Identification of 
misfolded SODl in sporadic ALS argues for a central role for 
SODl in ALS regardless of genetic or environmental initiators 
[19]. SODl'^'''''^ mice are an indispensable staple in the ALS field 
because they exhibit many etiological similarities to, and present 
with disease progression that has fidelity to, clinical ALS 
presentation. Zebrafish, with conserved CNS and motoneuron 
physiology/ genetics, have proven very useful to study the genetic 
relationships between fALS genes via loss-of-function analyses 
[20-27]. Furthermore, toxic gain-of-function zebrafish have been 
generated overexpressing mutant SODl and that recapitulate all 
hallmark attributes of fALS [28,29]. A substantial percentage of 
fALS remains to be explained, and genes that modify such 
susceptibility have potential to be influential in late-onset disorders 
with complex interactions of genetics and/or environmental 
factors. 

The proposed requirement for BMP/GDF6 family genes in 
vertebrate NMJs is supported by at least three previous lines of 
inquiry. First, disruption of intracellular trafficking led to NMJ 
deficits in zebrafish mutants of atlastin, and it was argued that this 
deficit was due to disrupted trafficking of BMP Receptors [30]. 
Second, BMPs and BMP Receptors are increased in expression at 
NMJs during recovery from experimentally induced traumatic 
injury in mice [5]. Third, a role for TGFfi/BMP/GDF signaling 
in ALS progression is suggested by the increase in pSMADs in 
ALS inclusions observed in both patient pathology and mouse 
models [31-33]. Although these lines of evidence are encouraging 
and suggestive, a role for BMP ligands in vertebrate neuromus- 
cular degeneration remains unaddressed. 

Our recent success in raising topical mutant zebrafish to 
adulthood [34] enabled us to address the questions above 
regarding a role for BMPs in the maintenance of vertebrate 
NMJs. Zebrafish harboring a nuU mutation in a homolog of GDF6 
were outcrossed, establishing a viable and fecund line of 
homozygous gdfSa mutants. Licreased cell death in at least some 
CNS neurons was reported during early development of these 
mutants [34]. This encouraged us to explore a potential role for 
gdf6a in neuromuscular degeneration, which would potentially 
implicate BMP/GDF6 genes as candidate loci in neuromuscular 
disease. We document progressive neuromuscular deficits in these 
mutants, including reduced muscle endurance, neuromuscular 
junction anomalies and loss of spinal motor neurons. Further, this 
mutation in gdf6a increased disease severity in a zebrafish model of 
ALS that over-expresses mutant SODl. These data integrate to 
compel GDF6 as a candidate locus that increases susceptibility to 
or progression of ALS. 



Results 

Gdf6~'~ zebrafish do not display motoneuron 
axonopathy in early development 

Because Gdf6 ' mice are not viable beyond embiyonic stages 
[34], we considered alternate animal models and asked if zebrafish 
deficient in a GDF6 homologue might present with neuromuscular 
deficits akin to ALS. Our recent work employed several 
generations of out-crosses to successfully isolate a viable line of 
homozygous null gdf6d^"'''"^ fish (ZFin ID: ZDB-ALT-050617- 
10) that can be raised to adulthood [34]. This nonsense mutation 
occurs in the pro-domain, at residue 55, abrogating production of 
mature GdfBa ligand and herein we refer to these mutants as 
gdf6a ' . These fish display microphthalmia [34], but are 
otherwise normal in development (Figure lA), adult morphology, 
and fecundity. A decreased longevity in gdf6a~'~ fish compared to 
their siblings was noted, not reaching statistical significance in a 
small cohort (Figure IB). 

To assess '\i gdf6a~'~ embryos display motoneuron disease, we 
utilized sensitive assays of motoneuron disease /ALS that are 
deployed frequently in embryonic zebrafish, consisting of exam- 
ining branching and pathfmding defects in GFP-positive primary 
motoneurons of transgenic fish [20-22,24-27]. gdf6a~'~ embryos 
did not display a significant increase in motoneuron pathfinding or 
axonopathies as determined in a transgenic background that 
enables sensitive detection of motoneuron morphology based on 
GFP fluorescence; both gdJ6a ' and siblings had a low rate of 
abnormalities in this scoring system (Supplemental Figure SI; 
p^O.315, nS:9 larvae per genotype). We confirmed that this assay 
was sensitive to ALS-related genetic lesions in our hands by 
delivering mRNA encoding human SODl with or without 
mutations associated with familial ALS (5.7 and 0.18 primary 
motoneurons affected per fish injected with SODl^*^ or 
SODl^^, respectively p<0.05), consistent with past results 
[26,27]. 

Considering ALS is a late-onset disease, it was of interest to 
characterize adult gdJ6a ' zebrafish with respect to ALS-like 
phenotypes. We examined fish at 9 to 18 months old (zebrafish are 
considered to be 'adult' at sexual maturity, ~3 months old and can 
thrive for 3-5 years of age) regarding muscle endurance, gait, 
motoneuron abundance and character of their neuromuscular 
junctions. 

gdf6a loss in adult zebrafish leads to neuromuscular 
disease, including deficits in endurance 

Deficits in muscle endurance is a hallmark of ALS progression 
consistently observed in murine and fish models [28]. We 
quantified muscle endurance in gdJ6a ' and sibling fish by 
assessing their ability to swim against a strong, accelerating water 
current [35], from which we determined the critical swimming 
speed (f/crit, the water flow velocity at which a fish can no longer 
maintain its position [35]; See Methods). At 9 months of age, 
gdfGa ' fish were found to have substantially reduced swimming 
endurance, with mean U„\i values approximately half those of 
gdf(^'* siblings (Figure 2A,B; p = 0.005, n=12 and 9 fish, 
respectively). At 18 months of age, the same significant difference 
was observed between genotypes (Figure 2B, p<0.0005, n = 4 and 
n = 3). Endurance was lower in the older fish within each 
genotype, though this difference between ages was not significantiy 
different when comparing within genotypes. 

The swim endurance (C/cjit) values in the wild type fish above are 
within the range of values obtained in past studies, being higher 
than [28], lower than [36,37] or comparable to [38] previous 
results. Differences between studies are likely due to methodology, 
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Figure 1. Zebrafish harboring homozygous mutations in gdfd can be raised to adulthood and do not display overt skeletal defects. 

A. gdf6^^^ fish are viable into adulthood, exhibit variably penetrant microphthalmia and normal body morphology. B. gdf6^'^ fish exhibit somewhat 
decreased survival compared to gdf6*^* siblings (n = 11 gdf6^^^ fish; n=4 gdf6*^* siblings). C, D. gdf6^^^ fish lack overt skeletal phenotypes, as 
revealed by (C) clearing and staining or by (D) microCT analysis. The latter is further represented as Supplemental Movie SI. Scale bars are 5 mm. A 
variety of fin morphologies were present in the fish examined, but these were neither different between experimental groups (genotypes) nor a 
significant covariant with swim performance (see Results). 
doi:10.1371/journal.pone.0089183.g001 



including duration of water velocity steps and flow chamber 
construction that can affect absolute values; thus comparisons 
between experimental configurations can be challenging, but the 
results are robust when compared within an experiment or 
laboratory. 

An alternative interpretation of our endurance data is that fish 
were less able to perform the .swimming test because of decreased 
availability of visual cues, due to microphthalmia. To rule out any 
confounding role for the visual deficits, swimming performance 
was assessed in a dark room under infrared light, which zebrafish 
are unable to detect (lighting conditions of the swimming 
endurance assay and zebrafish hght sensitivity are characterized 
in Supplemental Figure S2). A further alternative is that the fish 
have decreased endurance due to decreased activity patterns in 
their typical husbandry conditions, perhaps predicting decreased 
muscle fitness. Open field tests on groups offish were performed in 
conditions closely mimicking the husbandry conditions in which 
these fish were raised, and conducted in multiple replicates during 
each hour through a fuU circadian cycle. No significant difference 
in activity was found, arguing against mean activity as a 



confounding factor in these fish (Figure 2C, Supplemental Figure 
S2). The methods deployed were sensitive enough to detect lower 
average swim velocity in night vs. day (regardless of genotype, 
p<0.001), however no difference was detected between genotypes 
for average swim velocity during day (p = 0.247) or night 
(p = 0.814) (Figure 2C). 

A small difference was noted between genotypes in the open 
field test, only during mating activity at light onset (Supplemental 
Figure S2), which accounted for the small non-significant decrease 
in overall activity of the gdf6a~'~ fish. The difference observed 
during breeding is not relevant to the fish under consideration 
here, because they were raised and maintained in tanks that 
included sibling fish with normal eyes. When housed with 
normophthalmic siblings (our standard practice), gdf6a ' fish 
had ample activity during breeding times and coordinately robust 
fecundity; this contrasts groups comprised solely oi gdJGa ' fish 
that had no breeding success. Thus the open field test we deployed 
under-represented the gdf6a ' fish's activity during mating 
behaviour at light onset. In sum, the data argue against visual 
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Figure 2. Zebrafish harboring homozygous mutations in gdfS exhibit decreased endurance. A. gdf6~'~ fish have a lower endurance 
compared to gdf6*^* siblings as measured by increasing water velocity in a swim channel to determine the Uc„t i-e. 'critical swimming speed' which is 
taken to be the highest speed that a fish can swim at for a period of several minutes before exhaustion. Each fish tested at 9 months is plotted. B. The 
same data in panel A (9 month) plotted along with 1 8 month old siblings. gdf6^^^ fish have approximately 50% lower endurance compared to gdf6*^* 
siblings at each age (*p = 0.005, **p<0.01). Endurance trends towards being decreased in older fish of each genotype, but this difference is not 
significant. C. Open field tests of average swim velocity during 10 minutes of each hour through a circadian cycle, in tanks replicating lifetime 
husbandry conditions, show a lower mean movement at night (in either genotype, **p<0.01) but no difference between genotypes in any measure 
during day (p = 0.247) or night (p = 0.814). 
doi:1 0.1 371 /journal.pone.00891 83.g002 



system deficits playing either an acute or long-term role in the 
muscle deficits of microphthalmic gdf6a~^~ fish. 

A further alternative explanation for deficits in swimming 
performance could be morphological changes resulting from long- 
term loss of gdf6a. Skeletal changes might be anticipated from 
BMP namesake functions in bone morphogenesis, from skeletal 
defects observed in patients with GDF6 mutations [39-42] or past 
analysis of zebrafish gdf6a morphants [43]. In adult zebrafish, 
skeletal deficits in the axial skeleton were not observed (Figure 1, 
Supplemental Movie SI). Body morphology was also assessed, and 



no significant morphometric covariates of Z/crit were observed, 
including condition factor (Figure S2). Caudal fin length was not 
different between experimental groups (One-way ANOVA 
Fs,.')! = 1.320 p = 0.28) and was not a significant covariant in the 
relationship between experimental group and swimming perfor- 
mance (MLR: F] 3] =0.36 p = 0.55). Thus there was no evidence 
for a role of skeletal system malformation in the swimming deficits 
of gdf6a ' zebrafish. 
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gdf6a loss in adult zebrafish leads to neuromuscular 
disease, including disruption of neuromuscular junctions 

Histological hallmarks of ALS include neuromuscular deficits, 
comprising loss of motoneurons at the level of the spinal cord and 
abnormal structure of neuromuscular junctions (NMJ). Nine 
month old gdf6a~'~ zebrafish were noted to have disrupted 
NMJ, including increased mean volume and greater variation of 
the presynaptic motoneuron compartment compared to wild type 
siblings (Figure 3). This presynaptic motoneuron volume, when 
normalized to the post-synaptic volume, was nearly 3-fold greater 
in gdf6a ' mutants (p = 0.009, n = 5 fish per genotype). Co- 
localization of pre- and post-synaptic compartments was not 
significantly different based on genotype (Figure 3). Increased 
volume of pre-synaptic compared to post-synaptic NMJ compart- 
ments has also been observed in ALS model zebrafish overex- 
pressing mutant SODl [28]. 

Further, gdf6a ' fish at nine months of age were noted to have 
approximately 50% fewer spinal motoneurons compared to 
wildtype fish (p<0.001. Figure 3, n>50 sections from 4 fish per 
genotype), akin to changes observed in ALS model fish overex- 
pressing mutant SODl [28]. 

Disruption of gdf6 hastens disease progression in a 
zebrafish model of ALS 

We reasoned that \S gdf6a mutation hastens disease progression, 
or increases susceptibility to ALS, then this should be recognizable 
in ALS animal models lacking gdfSa function. We bred gdfSa 
mutants to an established zebrafish model of ALS: transgenic oslO 
zebrafish that over-express zebrafish Sodl*^'^"*^ and develop 
symptoms of ALS very reminiscent of transgenic mice over- 
expressmg human SODr'"*"^ [28]. Breeding gdf6a'~ ;Tg(.sodl:- 
sodlG93R,hsp70l:DsRed)o.sl() fish generated six genotypes, consist- 
ing of three gdf6a genotypes (+/+, +/— or —/—), each with the 
OS 10 transgene being present or absent. Herein we define 'bigenic' 
fish as those both homozygous nuU for gdf6a and possessing the 
oslO transgene, i.e. gdf6a ;Tg(sodl:sodlG93R,hsp70l:DsRed)oslO. 

Similar to the observations on fish larvae with gdf6a lesions 
alone, no significant increase in axonopathy was noted in primary 
motoneurons when oslO was also present in the larval genome 
(Supplemental Figure SI). In older fish, amongst the six genotypes 
examined, it was apparent that a subset of bigenic gdf6a~'~;osl0 
fish failed to thrive. The time to reach 75% survival was 
approximately 6 months in bigenic ALS model zebrafish lacking 
normal gdf6a function, but more than 13 months in all other 
genetic combinations, including the same ALS model fish with at 
least partial gdf6a function (Figure 4A). 

Muscle endurance was assessed in bigenic gdf6a ' IDS 10 fish 
and their siblings, using younger fish (4 months old) in an effort to 
capture any potential acceleration of disease progression. Deficits 
in swimming endurance attributable to gdf6a genotype alone were 
less dramatic at this young age: gdf6a~'~ fish (lacking the oslO 
transgene) had f/,.rit values approximately 75% those of their 
siblings but these differences were not statistically different 
(Figure 4B). This smaller endurance deficit in young fish contrasts 
the significant reductions observed at 9 and 18 months of age 
(where Uc,r\t was ~50% that of siblings. Figure 2B), providing 
additional support for the contention that gdJ6a loss of function 
leads to neuromuscular degeneration in a progressive manner. 

Similarly, young oslO fish consistently displayed decreased 
endurance within each gdJGa genotype, similar to what was 
observed previously for older oslO animals [28], though in no case 
did this rise to the level of statistical significance (Figure 4B). Only 
when the genetic lesions were combined were statistical differences 



attained, such that bigenic gdJ6a ;o.sl() fish had L'^jit values 
approximately 65% those of wUd type fish (Figure 4B, p<0.05 by 
MANOVA, sample size = 37, with 5-7 individuals per genotype as 
indicated on Figure abscissa). Thus, in young adult fish, gdf6a and 
oslO genotypes both imposed effects on muscle endurance that 
were consistent with neuromuscular disease, but these deficits only 
reached statistical significance when both genetic lesions were 
combined. 

Muscle strength/ power was also assessed in these six genotypes, 
by documenting the ability of the fish to sprint against increasing 
velocities of water flow. This constant acceleration test was similar 
to the Z7(:rit endurance test, but with water velocity increased at an 
accelerated rate such that fast-twitch white muscle became 
dominant, allowing the velocity at fatigue iU^.^ to be calculated 
[44,45] . Neither the presence of the oslO transgene nor loss oigdfGa 
function were individually potent in reducing swimming strength 
in these young fish (Figure 4C). Only when the two genetic lesions 
were combined was any statistically significant deficit in swimming 
power observed, compared to wild type fish. Thus bigenic gdf6a ' 
~ -fislO fish had f/max values approximately 75% those of wildtype 
fish (Figure 4C, p<0.05 by MANOVA, sample size = 37, with 5-7 
individuals per genotype as indicated on the Figure). 

These same bigenic gdf6a~ '~ ;oslO fish and their siblings were 
subsequently examined for neuromuscular junction abnormalities. 
Both oslO and bigenic gdf6a~^~ ;oslO fish at 7 months of age had 
increased presynaptic/postsynaptic volume ratios compared to 
wildtype siblings, though the differences did not reach statistical 
significance (Figure 5A, B. KruskaU-WaUis ANOVA p = 0.134, 
n = 6 WT, 4 oslO, 6 gdf6a~'~ ;oslO). Further, NMJs in bigenic 
gdf6a~^~ ;oslO had significantly decreased overlap of synaptic 
compartments compared to either oslO or wildtype siblings 
(Figure 5B) (KruskaU-Wallis ANOVA, p<0.001). Mander's 
colocalrzation coefficients for the postsynaptic compartment were 
significantiy lower for bigenic gdfSa ' ;oslO fish, implying loss of 
synaptic connectivity because postsynaptic compartments were less 
colocalized with presynaptic compartments despite the increase in 
presynaptic size (Figure 5B) (KruskaU-Wallis ANOVA, p<0.05). 
Thus the gdf6a ' mutation does not alter co-localization of NMJ 
pre- and post-synaptic compartments on its own (Figure 3C), but it 
further exacerbates the ALS-like NMJ abnormalities observed in 
bigenic fish expressing overexpressing mutant SOD 1 . 

Overall, loss of GdfB function led to an acceleration of disease 
progression in this ALS animal model, as measured by longevity, 
integrity of neuromuscuslar junctions, muscle endurance and 
muscle strength. These data are in good accordance with the 
hypothesis that deficits at the GDF6 locus can sensitize animals to 
onset and/ or progression of neuromuscular degeneration akin to 
ALS. 

Discussion 

Zebrafish harboring homozygous nuU mutations in a GDF6 
homologue were found to be viable, fecund and robust into 
adulthood, though exhibiting somewhat reduced longevity. 
Homozygous gdf6a ' fish had progressively severe deficits in 
endurance, loss of spinal motoneurons and disrupted neuromus- 
cular junctions compared to the wild type siblings, very much akin 
to ALS model zebrafish. Furthermore, progression of neuromus- 
cular disease in ALS model zebrafish was accelerated when this 
homologue of GDF6 was mutated: older zebrafish expressing 
mutant SODl exhibit ALS-like symptoms, but young fish only 
exhibited significant ALS-like neuromuscular deficits when on a 
gdJ6a~'~ background. 
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Figure 3. Zebrafish harboring homozygous mutations in gdf6 exhibit altered neuromuscular junctions and fewer spinal motor 
neurons. A-B. Assessment of neuromuscular junction (NMJ) morphology including presynaptic (synaptophysin) and postsynaptic (otBTX) 
compartments in 9 month old gdfS^^^ and gdf6*^* siblings by immunohistochemistry. C. ALS-like increases in motoneuron pre-synaptic volumes are 
observed in gdf6^'^ fish vj/hen normalized to post-synaptic volumes (**p = 0.009, n = 5 fish per genotype). Coefficients of colocalization for pre- and 
post-synaptic compartments are not altered, as is expected in later-stage zebrafish ALS models. D-E. IVIotor neuron cell bodies were identified in 
cross-sections of spinal cord using immunohistochemistry against choline acetyl transferase (ChAT, e.g. arrowheads) in nine month old gdf6*^* and 
gdf6^^^ fish (panels A and B, respectively). Bottom panels affirm motoneuron cell body identification using actin and nuclear counter-stains. 
F. gdf6~^~ fish have approximately 50% the abundance of spinal motor neurons compared to sibling gdfS*''* fish. (Mann-Whitney U Test, *** p<0.001, 
n>50 sections from 4 fish per genotype, researcher blinded to genotype during quantification). Scale bar = 60 |im in A,B and 50 jim in D,E. 
doi:1 0.1 371 /journal.pone.00891 83.g003 



Regarding ALS and other neuromuscular diseases, support for a 
BMPs' relevance is continuing to be revealed by genetic studies in 
Drosophila, including that a homologue of BMPs and homologues 
of BMP receptors are each required for normal NMJ development 
in flies [1,4,46,47], but suggestions of such functions for BMP 
pathway components rarely extend to vertebrates. Within 
vertebrates, BMP receptor type II has recently been locahzed to 
mouse NMJs, and one of its ligands BMP4 was found in close 



proximity [48]. BMP4, BMPS and their receptors have also been 
shown to have a role at some synapses in the CNS [49] . A role for 
BMPs in the vertebrate neuromuscular system has more frequently 
been suggested following motoneuron or spinal cord injury [50]. 
Recent work has also demonstrated a role for proper trafficking, 
including that of BMP receptors, in primary motor axon 
outgrowth during zebrafish development [30] , instigating a search 
for which component of the BMP gene family might underpin this 
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Figure 4. Mutations in gdfe sensitize SOD1^G93R zebrafisKi to 
develop ALS-lii<e symptoms. Six genotypes combining gdf6^^^ 
alleles and SOD1^G93R alleles were examined in endurance and sprint 
tests, which primarily measure red and white muscle respectively. The 
SOD1^G93R fish mimic many aspects of ALS progression [28]. 
Swimming performance was measured in all six genotypes (determined 
by outcrosses and examining progeny) at 4.5 months of age. A. bigenic 
fish, i.e. gdf6^^^ mutants expressing SOD1^G93R, had decreased 
survival compared to all other genotypes, i.e. their siblings { = gdf6*^* 
or gdfe^'^) with or without SOD1'^G93R. B. bigenic gdf6^'' mutants 



expressing SOD1^G93R had significantly lower endurance compared to 
heterozygous siblings also expressing SOD1'-G93R, and to WT and 
heterozygous siblings without SODT'GgSR. C. Sprint test demonstrates 
a significant deficit in white muscle function only when SOD1^G93R and 
gdf6^^^ genotypes are combined ( = 25% deficit). (ANOVA and multiple 
linear regression, p<0.05; sample size indicated below graph). 
doi:10.1371/journal.pone.0089183.g004 

effect. A potential example may be GDF15, a member of the BMP 
subfamily with divergent properties and typically associated with 
macrophage recruitment and inflammation, in that GDF15 
disruption causes loss of motoneurons in mice [51]. A link 
between BMP signaling and neuromuscular disease is suggested by 
the observation of phosphoSMAD proteins (downstream signaling 
components of TGFp/BMP/GDF signaling) being enriched in 
pathological ALS inclusions [31-33,52]. Human genetics suggests 
that loci associated with hereditary spastic paraplegia, including 
ATLASTIN, can regulate BMP receptor trafiicking and effect 
neuromuscular disease [30,50]. Despite this substantial level of 
interest, we are not aware of previously published data from 
vertebrate experiments or human genetics directly linking BMP/ 
GDF gene family members (ligands or receptors) with ALS or any 
neuromuscular disease. 

A substantial proportion of ALS susceptibility remains unex- 
plained with respect to human genetics. Analyses herein support 
GDF6, a member of the BMP family, as worthy of further 
consideration. A role for GDF6 homologues in progressive, late- 
onset neuromuscular disease was supported herein by muscle 
physiology, histopathology and behavioural assessment of zebra- 
fish carrying mutations in gdfSa. Further, gdf6a loss of function 
accelerated disease progression in a zebrafish model of ALS. The 
synthesis of this data combines with recent reports of a role for 
GDF6 in congenital and late-onset photoreceptor degenerations 
[34,53] to compel a role for GDF6 as a modifier gene in the 
etiology of disparate neuropathies. 

BMPs in neurodegenerative disease 

Inappropriate growth factor expression and signaling have been 
linked to inflammation and apoptosis in other neural diseases such 
as Alzheimer's, Parkinson's and, in one instance, Creutzfeldt- 
Jakob Disease, which are characterized by inappropriate accu- 
mulation of misfolded proteins, leading to neuron death via a 
multitude of proposed pathological pathways [54-56]. TGFfil and 
TGFP2 expression are high in brains afllicted with Alzheimer's, 
causing inflammation, astrogliosis and neuron death through such 
pathways as NADPH oxidase-induced oxidative stress and 
caspases [57,58]. Xenopus TGFfil, a member of the same gene 
superfamily as GDF6, is released from Schwann ceUs and can 
promote synaptogenesis at NMJ [59]. TGFfil is an effector of 
motoneuron survival (refs in [60]) and has been proposed as a 
mediator of increased ALS susceptibility in athletes abusing 
anabolic steroids [60]. Based on these similarities, it is plausible 
that GDF6 dysfunction could contribute to any of these pathways 
as they play out in ALS. Based on our data, we believe the 
expansion of ALS studies to the investigation of TGFfi signaling is 
warranted, beginning with GDF6. 

Therapeutic implications 

Application of BMPs as treatment has been investigated 
previously [50], especially foUowing injury, suggesting that further 
knowledge of GDF6 function could lead to strategies that replace 
BMP function. A number of TGFfi superfamily members are 
neuroprotective, such as BMP7 in cerebral ischemia and GDF5 in 
Parkinson's disease [61,62]. GDF 15 plays a role in the pathology 
of motoneuron loss, but is also required to effect the therapeutic 



PLOS ONE I www.plosone.org 



7 



February 2014 | Volume 9 | Issue 2 | e89183 



Gdf6 in Neuromuscular Degeneration 




B 



U1 

+1 
cu 

CD 
> 

■a 

M 

"to 

E 



1.5 



0.5 



□ wild type 

□ Sfc//6o*/^ sibling; SOD1<=93r 

□ gc//6o-/- mutant;SODl°93R 



* * * 




NMJ colocalization 
coefficient 



Ratio pre- /post- 
synaptic volume 



Figure 5. Disruption of gdf6 function exacerbates neuromuscular junction abnormalities in ALS model zebrafish. ALS model zebrafish 
possess disruptions to neuromuscular junctions (NMJ), and loss of gdf6 function exacerbates this by 7 months of age. A. The presynaptic junctions 
(labeled with synaptophysin antibody) and postsynaptic junctions (labeled with fluorescently tagged otBTX) in ALS model zebrafish expressing the 
mutant SODl'^'^'' show punctate morphology, deviations in presynaptic volume and less overall colocalization compared to WT sibling junctions. 
Some abnormalities are exacerbated in bigenic siblings expressing the mutant SODl'^'^'' that are also gdf6^^^ (scale bar is 40 \im). B. Quantification 
of these NIVIJs suggests the presynaptic/postsynaptic volume ratios of SODI^'^" and bigenic gdfd^^^ ■,SOD^'^^^'' zebrafish are larger than those of WT 
siblings at this age, though these differences do not rise to statistical significance (Kruskall-Wallis ANOVA, p = 0.134; n = 6,4,6 for WT, 5001°'^^ and 
bigenic fish respectively). Colocalization coefficients, that measure overall colocalization of presynaptic and postsynaptic junctions, are altered in 
these fish. The values for SGDI*^'^" zebrafish are significantly lower than wild type sibling values, indicating that presynapses and postsynapses 
overlap less, as characterized previously for this transgenic ALS model [28]. Bigenic SODl'^'^" zebrafish that are also gdf6-/- have a dramatically 
lower coefficient than either sets of siblings, including being 30% lower than ALS model SODI^^^" fish with normal gdf6a (*p<0.05; ***p<0.001. 
Kruskall-Wallis ANOVA with pairwise comparisons). 
doi:1 0.1 371/journal.pone.00891 83.g005 



effects of GDNF [51,63], illustrating how these signaling factors 
are able to participate in either neuroprotection or degeneration. 

Therapeutic potential also exists in small molecules such as 
aminopropyl carbazoles and their derivatives. For example it has 
recently been demonstrated that P7C3 can ameliorate phenotypes 
in zebrafish gdfSa mutants [34], and P7C3 has also recently been 
shown to be efficacious in mouse models of ALS and Parkinsons 
Disease [64,65]. Efficacy of P7C3 seems to be driven by a 
combination of pro-neural and neuroprotective signals, and 



identifying the mechanisms that underpin this could be aided by 
zebrafish models and/ or by investigation of its etiological effects in 
diverse neurodegeneratioiis. 

Finally, GDF6 may represent an entry-point into unexplored 
pathogenic mechanisms of vertebrate neuromuscular degenera- 
tion. In cases where such mechanisms have been identified, the 
etiology of ALS can be associated with a variety of cellular 
mechanisms. Ferraiuolo and colleagues [66] review these to 
include oxidative stress, RNA processing, endosomal trafficking. 



PLCS ONE I www.plosone.org 



8 



February 2014 | Volume 9 | Issue 2 | e89183 



Gdf6 in Neuromuscular Degeneration 



protein degradation, cytoskeletal integrity and glutamate excito- 
toxicity. It is conspicuous tliat, despite being a member of the BMP 
family with diverse functions, GDF6 is not currently associated 
with any of these etiological categories. 

Three mechanisms whereby BMP ligands and receptors affect 
NMJs can be suggested as potential routes to explore in the 
zebrafish system. First, retrograde signaling of BMP ligand from 
the muscle to the NMJ, and retrograde transport of receptors to 
activate transcription in the nucleus of motoneurons is important 
for motoneuron survival in Drosophila [4,50]. Second, BMP in 
Drosophila is required for stabilizing the NMJ through non- 
canonical signaling via LIM Kinase domain 1 [4,50]. Third, 
vertebrate BMP ligands have recendy been revealed to bind Agrin 
[67], and thus directly affect the clustering of acetylcholine 
receptors during formation and stabilization of the NMJ. In 
support of the latter we found that the NMJ morphology, as 
measured by a-bungarotoxin labelling of acetylcholine receptors, 
was disrupted in gdf6a ' mutants (Figures 3, 5), suggesting one 
priority for future analysis. 

Conclusion 

In summary, past conclusions regarding a role for Drosophila 
BMP ligands in effecting NMJ can be extended to vertebrates, 
supporting conservation of mechanisms that may be of substantial 
import to degenerative neuromuscular diseases. Zebrafish present 
a tractable system for assessing the role of BMP/GDF ligands 
(data herein) and their receptors [30] with respect to neuromus- 
cular disease. The data also synergize to argue for further analysis 
of GDF6 and other BMP genes as modifier loci, affecting 
susceptibility to ALS and perhaps a broader suite of neurodegen- 
erative diseases. Further efforts to assess the role of vertebrate 
BMPs in motoneuron disease are warranted. 

Methods 

Ethics statement 

All fish husbandry and procedures were completed with 
approval of the University of Alberta Animal Care and Use 
Committee: Biosciences, in accordance with the Canadian 
Council on Animal Care. 

Zebrafish husbandry and establishment of mutants and 
transgenic crosses 

Zebrafish bearing gdJGd'^'^^ ' (also known as radar '■'^''''■'-^^ 
[13], Zfin ID ZDB-ALT-050617-10), referred to here as gdf6~'~, 
previously identified as larval lethal, were recently acquired as 
viable adults through breeding of heterozygous mutants, preceded 
by backcrosses onto an AB background [34]. Homozygous 
gdJGd'^'^^''^^^ larvae were identified by the presence of microph- 
thalmia at 3dpf and were raised separately from their siblings. 
Mutants and siblings were confirmed in several cases by 
genotyping from finclip derived DNA at adulthood (using primers 
5'-TGTGAGACACGGCTCCACTT-3' and 5'-GCAGGACGA- 
GAGCTTACCAA-3') and/ or by examining the Mendelian 
frequency with which their progeny exhibited microphthalmia. 
At 9 and 18 months, GdfBa"^"^""^ adults and WT siblings 
underwent a swim channel assay and were subsequendy sacrificed 
for NMJ and spinal cord immunocytochemistry. 

ALS model zebrafish, expressing mutant SODl'"''* under 
control of the zebrafish sodl promoter [28], were generously 
provided by Christine E. Beattie. The fish, Tg(.sodl:- 
sodlG93R,hsp70l:DsRed)oslO (ZFin ID: ZDB-GENO- 10 1006-6), 
herein referred to as oslO, can be readily identified from their 
expression of DsRed under a promoter that is independent of the 



sodl expression cassette, but linked immediately contiguous to it in 
the construct. The osl 0 line was crossed to gdJGd'^^'"''^'^^ fish, and 
screened for the presence of osl 0 transgene following heat shock at 
2dpf [28]; the resultant compound heterozygous fish were inbred 
to acquire the combinations oi gdf6d'^'^ and oslO in the various 
combinations reported. 

Zebrafish husbandry and breeding used standard methods [68], 
with 14:10 L:D light cycle, 1300 ±100 |xS conductivity at 28±rC 
water temperature. Fish were fed a diet dominated by brine 
shrimp and trout chow. 

Assessing primary motoneuron axonopathy in 
embryonic zebrafish 

Zebrafish embryos were scored for abnormal primary moto- 
neuron axonopathy, as previously described [27], in transgenic 
animals where primary motoneurons abundantly express GFP 
under control of the promoter mnxl (also known as hb9, TStn 
ID: ZDB-GENO-080606-250). Breeding created gdJ6a^'~ ; 
Tg(mnxl:GFP)ml2 fish. These were bred to gdf6a' ;Tg(sodl:- 
sodlG93R,hsp70l:DsRed)oslO fish described above to generate 
genotypes including gdf6a fish and siblings with labelled 
motoneurons and/or expressing mutant Sodl''^''^. 

At 30hpf, larvae were visualized under a Leica MZ16F 
fluorescent stereomicroscope and axon branching was scored 
while the researcher was blind to genotype. As per established 
methods [27], branching of the GFP positive motoneurons was 
considered abnormal if it initiated at or above the ventral aspect of 
the notochord. Branching that was observed below this point was 
considered normal. The yolk sac extension was used as a 
landmark, and only the ten motoneurons immediately above it 
were scored for defects (Figure SI). Subsequent to quantifying its 
axonopathy, each of the larvae was reared individually in 24 well 
plates until they were old enough to assess genotypes based on 
microphthalmia, and heat-shocked for one hour prior to screening 
for DsRed expression associated with oslO. To assess the sensitivity 
of the axonopathy assay in our hands, Tg(mnxl:GFP)ml2 fish were 
injected with 900 pg of mRNA encoding wUdtype or A4V human 
SODl. niRNAs were synthesized and delivered as per previous 
methods [69]. Axonopathy scores per fish were analyzed using a 
Kruskall-WaUis ANOVA followed by a pairwise Mann-Whitney 
U-test, calculated in SYSTAT 12. All p-values reported are two- 
tailed and significance was set at p<0.05. 

Adult fish - assessing endurance in swim channel assay 
and activity in open field test 

Individual fish were tested for their critical sustained swimming 
(^:rit)) which was the maximum velocity at which they could swim 
for a set period in a 1 0 L modified Brett-style swim tunnel (Loligo 
Systems, DK). Fish were placed in the swim tunnel under blackout 
hoods with infrared lighting (see below for characterization of 
lighting) and allowed to acclimate for 30 minutes to a low velocity 
current (5 cm/s). The flow rate was then increased every 
10 minutes by 4 cm/s. Fish were monitored via infrared camera 
(SAV-CD120; Mateo, QE) and Elgato video capture software 
(Elgato Systems, CA), and were considered to have reached 
maximum velocity when they were no longer able to maintain 
position in the current and fell into a screen at the back of the 
swim chamber. After fish fatigued the test was halted and the fish 
removed following a cool-down swimming period of 1 0 minutes at 
low velocity. The four month old fish were allowed to recover in 
the chamber at 5 cm/ s for 45-minutes before completing a second 
test with shortened, one-minute steps to determine their maximum 
sprint swimming speed (finax)- The U„it and U„„ values were 



PLOS ONE I www.plosone.org 



9 



February 2014 | Volume 9 | Issue 2 | e89183 



Gdf6 in Neuromuscular Degeneration 



calculated using the Brett formula [70]. AH fish were measured for 
morphological traits that are potential covariates of swimming 
endurance [36,38]. Morphology measurements included: head 
length, muscle length, standard length, fork length, total length, 
body height, body width, caudal peduncle, and mass. Condition 
factor was determined using the standard formula from mass (g) 
and body length (BL; mm), calculated as lO'^-mass/BL''. 

Open field tests on zebrafish were carried out on groups of three 
fish per tank under conditions nearly identical to their normal 
husbandry conditions regarding temperature, water chemistry, 
water flow rates and lighting (except that IR lights were used to 
enable video recording in the dark). Fish were transferred from 
their housing tanks to test tanks of the same dimensions and were 
acclimated for one hour prior to the start of testing (N = 6 tanks 
with 3 fish each). Zebrafish motion was then recorded from above 
for 10-minutes every hour for 24-hours, using IR-lit cameras (as 
above) and video surveillance software fVisio-soft; Mateo, QE). 
Average swimming speeds were determined from video (15 
measures per second) for each fish using behavioural analysis 
software (Ethovision XT8.5; Noldus, NL). Infrared lighting was 
characterized using a calibrated spectrophotometer (USB4000- 
UV-VIS Ocean Optics, Dunedin FL) connected to a 455 micron 
optical fiber (QP450-2-XSR, Ocean Optics) held at the approx- 
imate position that fish maintained during trials. Spectra were 
recorded to Spectra Suite software (Ocean Optics). 

Differences in the day- and night-time swimming speed between 
gdf6 ' and wild-type fish were assessed using a two-way repeated 
measirres analyses of variance with a Holm-Sidak post-hoc test. 
Differences in Vi^-^.^ (BL/s) and U,rit (cm/s) between genotypes 
were analyzed using one-way ANOVA's with Holm-Sidak post- 
hoc tests. Holm-Sidak adjusted p-values were reported. AU 
swimming performance and open field test data were log 
transformed to achieve normality. AU data were tested for 
normality (Shapiro-Wilk) and equal variance (Levene's test) prior 
to the use of parametric statistics. Multiple linear regression 
analyses were also performed using standard body length and 
genotype as explanatory variables of swimming performance to 
ensure that the detected differences were independent of fish size. 
The statistical analysis of swimming performance and activity was 
conducted using SPSS 20.0 (SPSS IBM, Chicago, IL) and 
Sigmaplot 1 1 (Systat Software, CA). 

Adult zebrafish tissue preparation and 
immunocytochemistry 

To assess neuromuscular junctions (NMJ), muscles from adult 
zebrafish were removed, cut transversely into quarters (anterior to 
posterior) and frKcd in 4% paraformaldehyde/0.1 M PO4 for 
24 hours at 4°C. Samples were prepared for cryosectioning as per 
previous methods [71] by washing thrice in 5% sucrose/0.1 M 
PO4 (20 minutes), once in 12.5% sucrose/0.1 M PO4 (30 min- 
utes), and overnight in 20% sucrose/O.I M PO4. Samples were 
embedded and frozen in OCT media (Tissue-Tek 62550-12). 
16 |J,m sections were serially placed on Superfrost Plus microscope 
shdes (Fislierbrand 12550-15). AU sections used in immunocyto- 
chemistry were of the second-most anterior quarter of muscle. 

NMJs were visualized by labeling pre-synaptic terminals and 
post-synaptic terminals as reported previously [28]. Briefly, shdes 
were blocked for 60 minutes in 10% NGS/PBSTw and incubated 
in 1:50 rabbit anti-synaptophysin (Invitrogen 180130) in 2% 
NGS/PBSTw overnight at 4°C. Slides were washed with PBSTw 
and incubated in 1:500 AlexaFluor 488 chicken anti-rabbit 
(Invitrogen A21441)/l:100 AlexaFluor 555-tagged a-bungarotox- 
in (Invitrogen B35451)/2% NGS/PBSTw overnight at 4°C. 
Sections were viewed with a Zeiss LSM 700 confocal mounted 



on a Zeiss Axio Observer.Zl microscope and imaged with ZEN 
2010 (version 6.0) software (Carl Zeiss Microimaging). 

To quantify motoneuron abundance, spinal cords were 
dissected whole and fixed in 4% paraformaldehyde/0.1 M PO4 
for 24 hours at 4°C. The spinal cords were then prepared for 
cryosectioning as previously described. Transverse sections 
(20 |J,m) of each sample were serially placed on 6 Superfrost Plus 
microscope slides. Motor neurons were labeled as previously 
described [28]. Briefly, slides were incubated in 0.003% HjOj/ 
PBS for 20 minutes, blocked for 1 hour in 5% normal donkey 
serum (NDS)/1% DMSO/PBSTw, washed, and incubated in 
1:100 goat anti-ChAT (Chemicon AB 144P)/2% NDS/1% 
DMSO/PBSTw for 3 days at 4°C. Shdes were rinsed and washed 
in PBSTw, and incubated in 1:500 Alexafluor 555 or Alexafluor 
488 donkey anti-goat (Invitrogen A-21432 and A-11055)/2% 
NDS/1% DMSO/PBSTw overnight at 4°C. Actin was labelled 
using phaUoidin conjugated to AlexaFluor-488 (Irndtrogen 
#A12379). Nuclei were stained with TO-PRO-3 (Invitrogen, 
7^T3605). After three 30-minute washes with PBSTw, slides were 
coverslipped and imaged. Sections were viewed with a Zeiss LSM 
700 confocal mounted on a Zeiss Axio Observer.Zl microscope 
and imaged with ZEN 2010 (version 6.0) software (Carl Zeiss 
Microimaging). Powerpoint 2008 for Mac (Microsoft) was used to 
assemble figures following being merged and/or linearly manip- 
ulated for brightness and contrast in Photoshop CS3, Zen confocal 
software (Zeiss), and/ or Imaris x64 7.4.0 (Bitplane). 

Histological analysis 

NMJ volumes were measured using the voxel counter plugin for 
ImageJ 1.45 (Wayne Rasband, National Institutes of Health; 
http://rsbweb.nih.gov/ij/index.html) and colocalization was an- 
alyzed in Imaris x64 (version 7.4.0, Bitplane). Researcher was 
blinded to genotype prior to image analysis. Values were 
normalized to WT values and statistical analysis was performed 
with KruskaU- Wallis ANOVA (volume and colocalization mea- 
surements) and Mann-Whitney U test (colocalization measure- 
ments) on SYSTAT 12. 

Motoneurons in spinal cord cross-sections were quantified as 
per established methods [28] by averaging the number of 
motoneuron cell bodies, identified as ChAT-positive objects 
greater than 10 |a.m, per section. Researcher was blinded to 
genotype during quantification and analysis. Statistical analysis 
was performed with Kruskall-Wallis ANOVA in SYSTAT 12. 

Clearing and staining of skeletal elements 

Clearing and staining of cartilage and bone was performed 
using established protocols [72]. Briefly, adult fish were fixed 
whole in 4% paraformaldehyde overnight, rinsed in distilled water 
and preserved step-wise in 30%, 70% and 95% ethanol. The fish 
were then placed in alcian blue stain (cartilage) (Sigma-Aldrich 
A5268) in 30% acetic acid/70% EtOH for 6 hours and placed in 
saturated sodiiun borate overnight. Depigmentation was per- 
formed using 15% peroxide/85% KOH for 25-50 minutes or 
until pigment was removed. Fish were partially cleared in 30% 
saturated sodium borate with 1:250 tissue culture grade trypsin 
(VWR, CA97061-708)(1/16 teaspoon/30 niL). Bones were 
stained using a 1:1 000 dilution of stock alizarin red dye (Sigma- 
Aldrich A5533) in 1% KOH for 1 hour, or until sufficiently pink. 
Afterward, fish were fuUy cleared in trypsin solution and preserved 
step-wise in 30% and 70% glycerin/ 1% KOH. Fishes were finally 
stored and imaged in 100% glycerin. 
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Imaging skeletal elements in zebrafisli by microCT 

Intact fish were fixed in 4% paraformaldehyde/5% sucrose in 
0.1 M PO4 bufier for 24 hours and washed in 0.1 M PO4 pH 7.4. 
buffer 3-tinies for 20 minutes to remove fixative. Fish were 
incubated overnight in Lugol solution (Sigma-Aldrich; No. L- 
6146) with gentle agitation followed by three 20-muiute baths in 
0.1 M PO4 buffer. Tissue was dehydrated in a graded ethanol 
series and kept at — 20°C. To orient the specimen appropriately 
for scanning, it was partially embedded vertically in 2% agarose, 
caudal fin down, with the top half of a 1 5 mL Falcon tube serving 
as a mold. Imaging was performed on a SkyScanll74 (Broker; 
Konticli, BE) compact Micro CT (50 kV x-ray source and inbuilt 
1.3 MP cooled x-ray camera). Raw data obtained from the 
scanner was initially reconstructed using NRecon v. 1.6.6 (Bruker; 
Kontich, BE) and further processed to optimize viewing of 
structures of interest using Osirix v. 5.0.2 (Open-source DICOM 
Viewer). 360° rotating specimen videos were produced using 
Osirix V. 5.0.2 for viewing. 

Supporting Information 

Figure SI In larval zebrafish, mutations in gdfSa do not 
appreciably sensitize SODl'^G93R zebrafish to develop 
ALS-like symptoms. Four genotypes combining Gdf6 ' alleles 
and SOD1AG93R aUeles were examined in Tg(HB9:eGFP) 
zebrafish expressing GFP in the axons of primary motor neurons 
(PMN), or via immunohistochemistry. A. Bracket indicates 
position of axons quantified, magnified in B,C. B. Normal 
primary motor axons. C. An example of an abnormal PMN axon 
(arrow). D. Quantification of primary motor axons in 30 hour 
post-fertilization embryosshow no difference based on gdf6a 
genotype; Thus effects are not developmental, and accord with a 
late-onset phenotype. Results indicated no significant effect of 
Gdf6 on the presence of axonopathies (p&0.315, n&5 larvae per 
genotype), although the highest rate of PMN axon abnormalities 
were observed in gdf6a ' ;S0D1'^"^'^^'^ larvae. 
(TIF) 

Figure S2 Eliminating alternate hypotheses that might 
account for differences in swimming behaviour between 
gdf6a ' zebrafish and their siblings. Our overall conclu- 
sion is that gdJ6a mutants have deficits in endurance due to a 
progressive loss of spinal motor neurons and disrupted neuromus- 
cular junctions. Alternative explanations for these fish having 
reduced endurance are eliminated here. A. Body morphology 
was not significantly different based on genotype. 
Condition factor was determined using the standard formula from 
mass (g) and body length (BE, in mm), calculated as 
100000Xmass/BE'\ Condition factor did not vary based on 
genotype. Sample sizes (number offish) indicated at the bottom of 
graph. B. Acute deficits in vision cannot account for 
differences between microphthalmic gdf6a~'~ fish and 
their wildtype siblings. Infrared lighting conditions during 
behavioural tracking of zebrafish excludes a role for visual 
dysfunction in the assays of fish actiwty, power or endurance. 
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Movie SI gdf6~ '~ fish have normal skeletal morphology 
compared to wildtype siblings. See Fignre 1 regarding scale 
bars. 
(MOV) 

Acltnowledgments 

We thank Bradley Scott and Allison Murray for use of their MicroCT 
scanner and training. We received technical support from Neel Doshi and 
Aleah McCorry. Christine Beattie was generous in sharing her osl 0 line of 
transgenic zebrafish. We appreciate the generous constructive comments 
made by Simonetta Sipione and Kelvin Jones on an earlier draft of this 
manuscript. The Zebrafish International Resource Center (ZIRC) 
provided the transgenic line mnxligfp. 

Author Contributions 

Conceived and designed the experiments: MGD MJHG KBT WTA. 
Performed the experiments: MGD MJHG DEW TCZ KBT WTA. 
Analyzed the data: MGD MJHG DEW TCZ KBT WTA. Wrote the 
paper: MGD WTA. Critically edited the manuscript: MGD MJHG DEW 
TCZ KBT WTA. 



4. Bayat V, Jaiswal M, Bellen HJ (2011) The BMP signahng pathway at the 
Drosophila neuromuscular junction and its links to neurodegenerative diseases. 
Gurr Opm Neurobiol 21: 182-188. 

5. Ruschke K, Hiepen G, Becker J, Knaus P (2012) BMPs are mediators in tissue 
crosstalk of the regenerating musculoskeletal system. GcU Tissue Res 347: .^21— 
.544. 

6. Marques G (2005) Morphogens and synaptogenesis in Drosophila. J Nciu'obiol 
64: 417-434. 



PLOS ONE I www.plosone.org 



11 



February 2014 | Volume 9 | Issue 2 | e89183 



Gdf6 in Neuromuscular Degeneration 



7. Marques G, Haerry TE, Grotty ML, Xue M, Zhang B, et al. (2003) Retrograde 
Gbb signaling through the Bmp type 2 receptor wishful thinking regulates 
systemic P^RTa expression in Drosophila. Development 130: 5457—5470. 

8. Eaton BA, Davis GW (2005) LIM Kinase! controls synaptic stability 
downstream of the type II BMP receptor. Neuron 47: 695—708. 

9. Dion PA, Daoud H, Rouleau GA (2009) Genetics of motor neuron disorders: 
new insights into pathogenic mechanisms. Nat Rev Genet 10: 769-782. 

10. Nishimura AL, Mitne-Neto M, Silva HC, Richieri-Costa A, Middleton S, et al. 
(2004) A mutation in the vesicle-trafficking protein VAPB causes late-onset 
spinal muscular atrophy and amyotrophic lateral sclerosis. AniJ Hum Genet 75: 
822-831. 

11. Ratnaparkhi A, Lawless GM, Schweizer FE, Golshani P, Jackson GR (2008) A 
Drosophila model of ALS: human ALS-associated mutation in VAP33A suggests 
a dominant negative mechanism. Plos One 3: e2334. 

12. Graff JM (1997) Embryonic patterning: to BMP or not to BMP, that is the 
question. CeU 89: 171-174. 

13. Gosse NJ, Baier H (2009) An essential role for Radar (GdlBa) in inducing dorsal 
fate in the zebrafish retina. Proc Nad Acad Sci USA 106: 2236-2241. 

14. French CR, EricLson T, French DV, Pilgrim DB, Waskiewicz AJ (2009) Gdf6a is 
required for the ijiitiation of dorsal-ventral retinal patterning and lens 
development. Dev Biol 333: 37—47. 

15. Chang C, Hemmati-Brivanlou A (1999) Xenopus GDF6, a new antagonist of 
noggin and a partner of BMPs. Development 126: 3347—3357. 

16. Sasagawa S, Takabatake T, Takabatake Y, Muramatsu T, Takeshima K (2002) 
Axes establishment during eye morphogenesis in Xenopus by coordinate and 
antagonistic actions of BMP4, Shh, and RA. Genesis 33: 86-96. 

17. Rissi M, Wittbrodt J, Delot E, Naegeh M, Rosa FM (1995) Zebralish Radar: a 
new member of the TGF-beta superfamUy defines dorsal regions of the neural 
plate and the embryonic retina. Mech Dev 49: 223—234. 

18. Miyazono K, Maeda S, Imamura T (2005) BMP receptor signaling: 
transcriptional targets, regulation of signals, and signaling cross-talk. Cytokine 
Growth Factor Rev 16: 251-263. 

19. Pokrishevsky E, Grad LI, Yousefi M, Wang J, Mackenzie IR, et al. (2012) 
Aberrant localization of PTJS and TDP43 is associated with misfolding of SOD 1 
in amyotrophic lateral sclerosis. Plos One 7: e35050. 

20. Hewamadduma CA, Grierson AJ, Ma TP, Pan L, Moens CB, et al. (2013) 
Tardbpl splicing rescues motor neuron and axonal development in a mutant 
tardbp zebrafish. Hum Mol Genet 22: 2376-2386. 

21. Kabashi E, Bercier V, Lissouba A, Liao M, Brustein E, et al. (2011) FUS and 
TARDBP but not SODl interact in genetic models of amyotrophic lateral 
sclerosis. PLoS Genet 7: el002214. 

22. Laird AS, Van Hoecke A, Dc Muynck L, Timmcrs M, Van den Bosch L, et al. 
(2010) Progi'anulin is neurotrophic in vivo and protects against a mutant TDP- 
43 ijiduced axonopathy. Plos One 5: el 3368. 

23. Metz GA, Merkler D, Dictz V, Schwab ME, Fouad K (2000) Efficient testing of 
motor fianction in spinal cord injured rats. Brain Res 883: 165—177. 

24. Schmid B, Hruscha A, Hogl S, Banzhaf-Strathmann J, Strecker K, et al. (2013) 
Loss of ALS-associated TDP-43 in zebrafish causes muscle degeneration, 
vascular dysfunction, and reduced motor neuron axon outgrowth. Proc Natl 
Acad Sci U S A 1 10: 4986-4991. 

25. Vaccaro A, Patten SA, Ciura S, Maios C, Therrien M, et al. (2012) Methylene 
blue protects against TDP-43 and FUS neuronal toxicity in C. elegans and D. 
rerio. Plos One 7: e42117. 

26. Van Hoecke A, Schoonaert L, Lemmens R, Timmers M, Staats KA, et al. (2012) 
EPHA4 is a disease modifier of amyotrophic lateral sclerosis in animal models 
and in humans. Nat Med 18: 1418-1422. 

27. Lemmens R, Van Hoecke A, Hersmus N, Gcelen V, D'Hollander I, et al. (2007) 
Overexpression of mutant superoxide dismutase 1 causes a motor axonopathy in 
the zebrafish. Hum Mol Genet 16: 2359-2365. 

28. Ramesh T, Lyon AN, Pineda RH, Wang C Janssen PM, et al. (2010) A genetic 
model of amyotrophic lateral sclerosis in zebrafish displays phenotypic hallmarks 
of motoneuron disease. Dis Model Mech 3: 652—662. 

29. McGown A, McDearmid JR, Panagiotaki N, Tong H, Al Mashhadi S, et al. 
(2013) Early interneuron dysfunction in ALS: insights from a mutant sodl 
zebrafish model. Ann Neurol 73: 246-258. 

30. Fassier C, Hutt JA, Scholpp S, Lumsden A, Giros B, et al. (2010) Zebrafish 
atlastin controls motility and spinal motor axon architecture via inhibition of the 
BMP pathway. Nat Neurosci 13: 1380-1387. 

31. Nakamura M, Ito H, Wate R, Nakano S, Hirano A, et al. (2008) Phosphorylated 
Smad2/3 immuno re activity in sporadic and familial amyotrophic lateral 
sclerosis and its mouse model. Acta Neuropathol 115: 327—334. 

32. Nakamura M, Kaneko S, Ito H, Jiang S, Fujita K, et al. (2013) Activation of 
transforming growth factor-beta/ Smad signaling reduces aggregate f ormation of 
mislocalized TAR DNA-binding protein-43. Neurodegener Dis 11: 182—193. 

33. Nakamura M, Kaneko S, Wate R, Asayama S, Nakamura Y, et al. (2012) 
Regionally different immunoreactivity for Smurf2 and pSmad2/3 in TDP-43- 
positive inclusions of amyotrophic lateral sclerosis. Neuropathol Appl Neurobiol. 

34. Asai-Coakwell M, March L, Dai XH, Duval M, Lopez I, et al. (2013) 
Contribution of growth differentiation factor 6-dependent cell survival to early- 
onset retinal dystrophies. Hum Mol Genet 22: 1432-1442- 

35. Tierney KB (2011) Swimming performance assessment in fishes. J Vis Exp. 

36. Plaut I (2000) Effects of fin size on swimming performance, swimming behaviour 
and routine activity of zebrafish Danio rerio. J Exp Biol 203: 813—820. 



37. McClelland GB, Craig PM, Dhekney K, Dipardo S (2006) Temperature- and 
exercise-induced gene expression and metabolic enzyme changes in skeletal 
muscle of adult zebrafish (Danio rerio). J Physiol 577: 739-751. 

38. Gilbert MJ, Zerulla TC, Tierney KB (2014) Zebrafish (Danio rerio) as a model 
for the study of aging and exercise: Physical ability and trainability decreases 
with age. Exp Gerontol50C: 106-1 13. doi: 10.1016/j.exger.2013.1 1.013. 

39. Asai-Coakwell M, French CR, Ye M, Garcha K, Bigot K, et al. (2009) 
Incomplete penetrance and phenotypic variability characterize GdfB-attribut- 
able oculo-skeletal phenotypes. Hum Mol Genet 18: 1110—1121. 

40. Asai-CoakweU M, French CR, Berry KM, Ye M, Koss R, et al. (2007) GDF6, a 
novel locus for a spectrum of ocular developmental anomaUes. Am J Hum Genet 
80: 306-315. 

41. Tassabehji M, Fang ZM, HUton EN, McGaughran J, Zhao Z, et al. (2008) 
Mutations in GDF6 are associated with vertebral segmentation defects in 
Klippel-Feil syndrome. Hum Mutat 29: 1017-1027. 

42. Setde SH, Jr., Rountree RB, Sinha A, Thacker A, Higgins K, et al. (2003) 
Multiple joint and skeletal patterning defects caused by single and double 
mutations in the mouse GdfB and Gdf5 genes. Dev Biol 254: 116-130. 

43. Ye M, Berry-Wynne KM, Asai-Coakwell M, Sundaresan P, Footz T, et al. 
(2010) Mutation of the bone morphogenetic protein GDF3 causes ocular and 
skeletal anomalies. Hum Mol Genet 19: 287-298. 

44. FaiTcU AP (2007) Comparisons of swimming performance in rainbow trout using 
constant acceleration and critical swimming speed tests. Journal of Fish Biology 
72: 693-710. 

45. Tierney KB (2011) Behavioural assessments of neurotoxic effects and 
neurodegeneration in zebrafish. Biochim Biophys Acta 1812: 381—389. 

46. Aberle H, Haghighi AP, Fetter RD, McCabc BD, Magalhaes TR, et al. (2002) 
wishful thinking encodes a BMP type II receptor that regulates synaptic growth 
m Drosophila. Neuron 33: 545-558. 

47. James RE, Broihier HT (2011) Crimpy inhibits the BMP homolog Gbb in 
motoneurons to enable proper growth control at the Drosophila neuromuscular 
junction. Development 138: 3273-3286. 

48. Chou HJ, Lai DM, Huang CW, McLennan IS, Wang HD, et al. (2013) BMP4 is 
a peripherally-derived factor for motor neurons and attenuates glutamate- 
mduced excitotoxicity in vitro. Plos One 8: e5844L 

49. Xiao L, Michalski N, Kronander E, Gjoni E, Genoud C, et al. (2013) BMP 
signaling specifies the development of a large and fast CNS synapse. Nat 
Neurosci 16: 856-864. 

50. HenriquezJP, Ki'ull CE, Osses N (201 1) The Wnt and BMP families of signaling 
morphogens at the vertebrate neuromuscular junction. Int J Mol Sci 12: 8924— 
8946. 

51. Strelau J, Strzclczyk A, Rusu P, Bcndiier G, Wiese S, et al. (2009) Progressive 
postnatal motoneuron loss in mice lacking GDF-15. J Neurosci 29: 13640— 
13648. 

52. Katsuno M, Adachi H, Banno H, Suzuki K, Tanaka F, et al. (2011) 
Transforming growth factor-beta signaling in motor neuron diseases. Curr 
Mol Med 11: 48-56. 

53. Zhang L, Lim SL, Du H, Zhang M, Kozak I, et al. (2012) High temperature 
requirement factor Al (HTRAl) gene regulates angiogenesis through trans- 
forming growth factor-beta family member growth differentiation factor 6. J Biol 
Chem 287: 1520-1526. 

54. Vawter MP, DiUon-Carter O, Tourtellotte WW, Carvcy P, Freed WJ (1996) 
TGFbetal and TGFbeta2 concentrations are elevated in Parkinson's disease in 
ventricular cerebrospinal fluid. Exp Neurol 142: 313-322. 

55. Li C, Zhao R, Gao K, Wei Z, Yin MY, et al. (201 1) Astrocytes: miplications for 
neuroinflammatory pathogenesis of Alzheimer's disease. Curr Alzheimer Res 8: 
67-80. 

56. Deininger M, Mcyermann R, Schluesener H (1995) Detection of two 
transforming growth factor-beta-related morphogens, bone morphogenetic 
proteins-4 and -5, in RNA of multiple sclerosis and Crcutzfeldt-Jakob disease 
lesions. Acta Neuropathol 90: 76—79. 

57. Wyss-Coray T, Lin C, von Euw D, Masliah E, Muckc L, et al. (2000) 
Alzheimer's disease-like cerebrovascular pathology in transforming growth 
factor-beta 1 transgenic mice and functional metabolic correlates. Ann N Y Acad 
Sci 903: 317-323. 

58. Hashimoto Y, Chiba T, Yamada M, Nawa M, Kanekura K, ct al. (2005) 
Transforming growth factor beta2 is a neuronal death-inducing ligand for 
amyloid-beta precursor protein. Mol Cell Biol 25: 9304—9317. 

59. Feng Z, Ko CP (2008) Schwann cells promote synaptogenesis at the 
neuromuscular junction via transforming growth factor-betal. J Neurosci 28: 
9599-9609. 

60. Galbiati M, Oncsto E, Zito A, Crippa V, Rusmini P, et al. (2012) The anabolic/ 
androgenic steroid nandrolone exacerbates gene expression modifications 
induced by mutant SODl in muscles of mice models of amyotrophic lateral 
sclerosis. Pharmacol Res 65: 221-230. 

6 1 . Perides G, Jensen FE, Edgecomb P, Rueger DC, Charness ME ( 1 995) 
Neuroprotective effect of human osteogenic protein- 1 in a rat model of cerebral 
hypoxia/ischemia. Neurosci Lett 187: 21—24. 

62. Hurley FM, Costello DJ, Sullivan AM (2004) Neuroprotective effects of delayed 
administration of growth/diflerentiation factor-5 in the partial lesion model of 
Parkinson's disease. Exp Neurol 185: 281-289. 

63. Krieglstein K, Strelau J, Schober A, Sullivan A, Unsicker K (2002) TGF-beta 
and the regulation of neuron survival and death. J Physiol Paris 96: 25-30. 



PLOS ONE I www.plosone.org 



12 



February 2014 | Volume 9 | Issue 2 | e89183 



Gdf6 in Neuromuscular Degeneration 



64. Tesla R, Wolf HP, Xu P, Drawbridge J, Estill SJ, ct al. (2012) Neuroprotective 
efficacy of aminopropyl carbazoles in a mouse model of amyotrophic lateral 
sclerosis. Proe Nad Acad Sci U S A 109: 17016-17021. 

65. De Jesus-Cortes H, Xu P, Drawbridge J, Estill SJ, Huntington P, ct al. (2012) 
Neuroprotective efficacy of aminopropyl carbazoles in a mouse model of 
Parkinson disease. Proe Nad Acad Sci U S A 109: 17010-1701."). 

66. Ferraiuolo L, Kirby J, Grierson AJ, Sendtner M, Shaw PJ (2011) Molecular 
pathways of motor neuron injury in amyotrophic lateral sclerosis. Nat Rev 
Neurol 7: 616-630. 

67. Banyai L, Sondereggcr P, Patthy L (2010) Agrin binds BMP2, BMP4 and 
TGFbetal. Plos One 5: cl0758. 

68. Westerlield M (2000) The zebrafish book. A guide for the laboratory use of 
zebralish (Danio rerio). Eugene, OR: University of Oregon Press. 



69. Kaiser DM, Acharya M, Leighton PL, Wang H, Daude N, et al. (2012) Amyloid 
beta precursor protein and prion protein have a conserved interaction aftecting 
cell adlresion and CNS development. Plos One 7: e51305. 

70. Brett JR (1964) The respiratory metabolism and .swimming performance of 
young sockeye salmon. Journal of the Fisheries Research Board of Canada 2 1 : 
1183-1226. 

71. Eraser B, Duval MG, Wang H, Allison WT (2013) Regeneration of cone 
photoreceptors when cell ablation is primarily restricted to a particular cone 
subtype. Plos One 8: e.')5410. 

72. Taylor WR, van Dyke GC (1985) Revised procedures for staining and clearing 
small fishes and other vertebrates for bone and cartilage study. Cybium 9: 107— 
119. 

73. AUison WT, Hamibcrgcr TJ, Hawryshyn CW, Temple SE (2004) Visual 
pigment composition in zebralish: Evidence lor a rhodopsin-porphyropsin 
interchange system. Vis Neurosci 21: 945-952. 



PLOS ONE I www.plosone.org 



13 



February 2014 | Volume 9 | Issue 2 | e89183 



